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Abstract 

The dominantly orbital state method allows a semiclassical description of quantum systems 
and is expected to be valid for large values of the orbital angular momentum and small radial 
excitations. At the origin, it was developed for two-body relativistic systems. Here, the method is 
extended to treat two- and three-body Hamiltonians in D > 2 dimensions with arbitrary kinetic 
energy and potential. It is shown that this method can produce quite accurate results, in some 
cases for the whole spectrum. 
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I. INTRODUCTION 



The dominantly orbital state (DOS) method is a technique to compute approximate 
solutions for quantum eigenvalue problems The first step is to compute the energy 
associated with a classical circular motion and then quantize the orbital angular momentum. 
This kind of semiclassical approximation is expected to be valid only for high values of this 
quantum number. Then, the radial motion is quantized around the physical circular orbit. 
This implies that the radial excitation must be small compared to the orbital excitation. As 
it will be shown, this method is very simple to implement. The only mathematical difficulty 
is the computation of the physical radius which is the solution of a transcendental equation. 
In favorable cases, the solution can be analytical. 

At the origin, the method has been developed for two-body relativistic Hamiltonians [l| 
and has been used to study hadronic systems l|-[3j]. The purpose of this work is to generalize 
this technique to treat two- and three-body Hamiltonians in D > 2 spatial dimensions with 
arbitrary kinetic energy and potential. This is motivated by the existence of non-standard 
kinetic energies in some physical problems, for instance in atomic physics with non-parabolic 



dispersion relation (4 1 , in hadronic physics with particle masses depending on the relative 
momentum |5j, or in quantum mechanics with a minimal length jf], Q]. Moreover, problems 
in D dimensions can appear in various physical situations. In particular, D = 2 systems can 

n 

be used as toy models for D = 3 systems [8j or are the natural framework for the physics 
of anyons js]. Three-body calculations are much more involved than two-body ones. The 
purpose of the extension of the DOS method to this kind of calculations is not to compete 
with numerical methods such that expansions in harmonic oscillator basis [10] or in Gaussian 



states 



but to provide rapidly approximate but reliable, and even analytical, results. 



This paper is organized as follows. Some elements for the physics in D dimensions are 
briefly reminded in Sec. [Til The DOS method for two-body systems is generalized in Sec. IIHI 
to .D-dimensional cases, and it is developed for three-body systems in Sec. IIV1 In both 
sections, applications are given to test the validity of the method. Concluding remarks are 
given in Sec. [V] 
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II. KINEMATICS IN D DIMENSIONS 



Let p and r be conjugate variables, then the operator p 2 in D > 2 dimensions can be 



written as 



12| 



2 2 i 2 f d 2 D — 1 d , 

P =Pr + ^ Where = - ( IJi + — — ) , (1) 



\ dr 2 r dr 



and where the orbital operator A 2 is such that [lJ 

A 2 Y lt{li} (Sl D ) = 1(1 + D - 2) Y h{pi} (n D ). (2) 

Functions ^/^(^d) are hyperspherical harmonics on the D-sphere with spherical coordi- 
nates Qjj. I is the orbital quantum number and {/i} are the magnetic quantum numbers. 
In semiclassical calculations, it is usual to replace the orbital factor 1(1 + a) by (I + b) 2 , 
the two forms being equivalent for / ^> 1, that is to say a domain of values for which the 
semiclassical approximation is expected to be relevant. In this case, we have 

A 2 D — 2 

P 2 =P 2 r + ^T With A = / + — — . (3) 

r 1 2 

III. TWO-BODY SYSTEMS 
A. General formula 

Let us consider the two-body Hamiltonian (H — c—1) 

H = T(p) + V(r), (4) 

with p = \p\ and r = \r\. For a circular motion (p r = 0), the semiclassical (A is quantized) 
energy is given by 

E(r)=T(^+V(r). (5) 

As usual, we assume that the physical orbit is determined by a minimum energy condition. 
This occurs for the radius tq given by 

-T'(-) =r V'(r ). (6) 
V r o/ 



potential 



14 




in 


16 



remarked in 16]. The corresponding energy is E = E(tq). 

3 



In order to compute the radial excitation, we will quantize the radial motion Ar around a 
circular orbit with a fixed radius r = r — Ar. This motion is controlled by the Hamiltonian 



AH = T[ x jpl + (ro+ A2 Ar)2 j + V(r + Ar) - E (7) 

For small values of the radial quantum number n, we can expand AH in powers of p 2 . and 
Ar. Keeping only the first non vanishing contributions, we obtain 

AH "T-/r + -2^ (8) 



" = ^7iy md k ^ ( 2roT '(^) +XT "^ ]] + v " {ro) - (9) 



with 



ro / 

The term in Ar is canceled by the circular condition for the orbit. Since / 3> 1, the 
radius r is expected to be large. Under this condition, we can assume that p 2 . ~ —d 2 /dr 2 , 
and AH reduces to a one-dimensional harmonic oscillator whose eigenvalues are given by 
AE = * (n+ |). Finally, the total energy, E = E + AE, can be written 



with -jT' ( — J = V\r ) and A = I + ■ (10) 



This system is expected to be valid when I S> 1 and n<l. The relevance of this formulation 
is tested in Sec. IHICI for three analytical solutions of f JTOj) . 



B. Case D = 2 and / = 



In the particular case D = 2 and I = 0, A = and the system (ITU]) is ill-defined. In 
principle, this is not a problem since the approximation is only valid for I ^> 1. Nevertheless, 
for some systems, it is possible to obtain a good accuracy even for small values of I (see 
following section). If the determination of ro is analytical, a generic value of A can be kept 
through the calculations. In the final formula, one must then check that the limit A — > is 
relevant. If only numerical calculations are possible to fix r , it is then interesting to dispose 
of an alternative method to treat specifically the case A = 0. 

4 



The WKB method is a powerful one to solve eigenequations without centrifugal term [17] . 
With our notation, the eigenvalue E can be computed by solving the following equation 

2 T- X (E-V(r)) dr = an + b, (11) 
Jo 

where a and b are two constants. With a vanishing orbital momentum for D = 2 (A = 0), 
the classical motion runs from to the turning point r* given by 

n = V- 1 (E-T(0)). (12) 

We can fix a and b in order that the WKB solution gives the exact result for the harmonic 
oscillator. We have then 

rV- l (E-T(0)) / i\ 

J T- 1 (E - V{r)) dr = n\n + -\. (13) 

Equations ffTUj) and f[T3"j) allow a determination of the energy spectrum for any value of A. 



C. Applications 

For the D-dimensional nonrelativistic harmonic oscillator, T(p) = p 2 / (2 m) and V(r) = 
kr 2 /2, the system fTTOjl gives 

E=J^- (2n + l + ^), (14) 



m \ 2 / 

which is the exact result 12J. By construction, this result is also obtained by ffT3|) for D = 2 
and I = 0. 

For the .D- dimensional Coulomb problem, T(p) = p 2 /(2m) and V(r) = —a/r, the sys- 
tem (ITUj) gives 

E = — ^- 2 + ma2 as) 



which are the first terms of the expansion of the exact solution 12] 



m a 2 



-^exact — — — ; D-l\ 2 



2(n + l + ^) 

for I 3> n. In the case D = 2 and 1 = 0, fTT3|) gives the exact result ([16 
Let us now look at the following semirelativistic Hamiltonian 



H = 2Vp2 + ar, (17) 



commonly used for the study of mesons composed of light quarks [18j]. The eigenvalues 
of (|17p are the masses of the system. An analytical solution of the system ffTOl) can also 
be obtained in this case. As it is remarked in [l| for the case D = 3, a more accurate 
approximation is obtained by computing the square energy E 2 but by dropping the term in 
(n + |) 2 , which is coherent with the DOS approximation. The final result is then 

E 2 = 8 a \V2 n + I + - ~ ^ + ^ j . (18) 

The Hamiltonian (ITT)) leads to Regge trajectories: The square mass is a linear function of I 
and n with slopes independent of D. The ratio between the radial and orbital slopes is \/2, 



which is in agreement with the value 7r/2 found in 19J for D = 3 with a Bohr-Sommerfeld 



quantization procedure. For D = 2 and 1 = 0, (TT3)) gives 

E 2 = Aira(n + ^ , (19) 

which is in quite good agreement with (ITS)) . The accuracy of formula (ITSl is tested in table fl] 
for D = 3, where the DOS solutions are compared with the eigenvalues computed with the 
high accuracy numerical Lagrange-mesh method [20J. One can see that the agreement is 
very good, even for small values of I and large values of n. It is expected that the method 
be more accurate for Hamiltonians close to the harmonic oscillator one, whose spectrum is 
exactly reproduced. 

IV. THREE-BODY SYSTEMS 
A. General formula 

Let us now consider a classical Hamiltonian for three identical particles interacting via 
the one-body U and two-body V interactions 

3 3 3 

H = J2T(\p t \) + J2u(\r l -R\)+ VQn-rjl), (20) 

i=l i=l i<i=l 

where Y^=i Pi = ® an d R — \ Yl^=i r i ^ s ^ ne center of mass position. A set of internal 
coordinates (x = X\ = r\ — r 2 ,x 2 = r 2 — r 3 ,x 3 = R) can be defined together with its 
conjugate variables (q = qi,q2,q?,)- In order to apply the DOS approximation, we will 



TABLE I. Eigenmasses Ej y/a of the Hamiltonian f)17|) for D = 3. For each set of quantum numbers 
{/, n}, the first line is an accurate value obtained by the Lagrange-mesh method [20], and the second 
line is given by (fT8|) . 



n = 1 



I = 3.157 4.709 5.889 6.871 

3.108 4.579 5.682 6.603 

1 4.225 5.457 6.483 7.375 
4.202 5.382 6.347 7.183 

2 5.079 6.130 7.047 7.867 
5.065 6.080 6.949 7.720 

3 5.811 6.724 7.577 8.338 
5.801 6.706 7.502 8.222 



impose a severe constraint on the system: the three particles move on a circular motion 
around the center of mass at the apex of an equilateral triangle. Under these conditions, we 
have x = \x\ = |r» — rj\ = v3|rj — R\ and q = \q\ = \pi\, V i and j. The Hamiltonian ( 12 Op 
can then be rewritten 

#dos = 3 [T{q) + W{x)\ with W{x) = U ( ^= j + V(x). (21) 



In this Hamiltonian, q 2 = + 

In order to interpret correctly the quantum numbers associated with the conjugate vari- 
ables x and q, once quantized, let us look at the total orbital angular momentum. In the 
center of mass frame, it is given by L = X^=i( r « — R) x Pi, and is such that 

\L\ = V3\xxq\. (22) 

If A (A) is the orbital angular momentum associated with L (x x q), we have A = 
Physically, A must be associated with the contributions brought by the two internal variables, 
and we can expect that A = (li + ^y^) + (h + ^j^) , where li and l 2 are the orbital angular 
momentums associated with these variables. Finally, we identify 

A = (L + D - 2) with L = h + l 2 . (23) 
V3 



In a semiclassical approach, the quantification of the radial motion for the three-body system 
is given by the integral of the quantity 3p r Ar, where p r is the radial momentum for one 
particle and Ar is the variation of the radius of the circular motion. We have 3 p r Ar = 
y3q r Ax, where q r and Ax are the equivalent quantities for the variables q and x. If 
v (n + 1/2) is the radial quantum factor associated with 3p r Ar (q r Ax), we have v = 
\/3(n + 1/2). Physically, v must be associated with the contributions brought by the two 
internal variables, and we can expect that v = \n\ + §) + (ri2 + |) , where ri\ and n 2 are the 
radial quantum numbers associated with these variables. Finally, we identify 

n + -= — ~L~*" with N = n 1 + n 2 . (24) 
2 v3 

Using the same procedure as in Sec. IIII A] we can find the DOS solution for f[2~T]) . It is 
written 



3 T (x ) + W ^ + y % l T (x ) + xl T (x ) T (x ) + A T (x ) W ^ v 7 ^ 

with -4 T ' ( — ) = W'(x ) and \ = -^(L + D-2). (25) 
x o \ x oJ v3 

Again, this system is expected to be valid when L ^> 1 and N ^ L. The relevance of this 

formulation is tested in Sec. HVCI for two analytical solutions of ( 121)]) . 

If the masses of the three bodies are different, the assumption \q\ = \pi\ is no longer 

justified. So, the present approximation scheme should not be applied. However, if a strong 

mass asymmetry exists between the particles, the problem can be reduced to a two- or 

one-body problem thanks to an adiabatic approximation. 



B. Case D = 2 and L = 

Again, if D = 2 and L = 0, the system ( 125]) is no longer well-defined. Using, the same 
considerations as in Sec. IIII B] the WKB method can be used to treat this particular case. 
This finally gives 

W/-i(f-T(0)) (E X 

l T-^--W(x)j dz = -j=(N + l). (26) 



From 



2l[ |. it is easy to determine that the exact solution for the general nonrelativistic D- 



dimensional three-body harmonic oscillator, T(q) = q 2 /(2m), U(x) = kx 2 and V(x) = px 2 , 



is written 



It can be checked that fl26l) gives the exact result with D — 2, li — l 2 — 0, and N = rii + n 2 . 



E= J — (k + 3p)(2n 1 + 2n 2 + l 1 +l 2 + D). (27) 
m 



C. Applications 

For the general three-body harmonic oscillator defined in the previous section, the sys- 
tem (l2"o"|) gives 

E = J-{k + 3p){2N + L + D). (28) 
V m 

which is the exact result with the identification N = ri\ + n 2 and L = l\ + l 2 . This confirms 
the interpretation chosen for the quantum numbers N and L in Sec. IIV Al By construction, 
this result is also obtained by f )26|) for D = 2 and L = 0. 

Let us now look at the following semirelativistic Hamiltonian 

3 3 3 

i=l i=l *<J=1 

used for the study of baryons composed of light quarks [22|. It can be seen as the three- 
body generalization of (I17p . The eigenvalues of (1291) are the masses of the system. An 
analytical solution of the system (|25p can also be obtained in this case. As for the two-body 
Hamiltonian, a more accurate approximation is obtained by computing the square energy 
E 2 but by dropping the term in (N + l) 2 . The final result is then 

E 2 = 12c(V2N + L + D -2 + v 7 ^) with c = a + V3b. (30) 

By comparing this equation with (iTSj) . one can see that the same ratio between the Regge 



slopes is predicted for two- and three-body systems, as expected 
(ESD gives 



22| . For D = 2 and L = 0, 



E = 6nc(N + 1), (31) 



which is in quite good agreement with (130]) . 

The eigenstates of the Hamiltonian ( 129|) have also been computed with a high accuracy 



method relying on the expansion of trial states on a harmonic oscillator basis [10] . The lowest 
eigenstates are characterized by a large component (generally more than 90%) spanned by 
one or several harmonic oscillator basis states with the same fixed value of iV = n\ + n 2 



and L = l\ + li- So they can be labeled by the number of quanta Q = L + 2 N, and the 
corresponding eigenvalues can be easily compared with the predictions of the DOS method. 
All eigenvalues for the completely symmetrical and the completely antisymmetrical states 
are presented in Fig. [TJ For the clearness of the figure, only the eigenvalue corresponding 
to the mixed symmetry state for Q = 1 is shown. One can see that the agreement is quite 
good, even if some degeneracies remain. In a pure harmonic oscillator picture, all masses 
associated with the same value of Q should be degenerate. It can be observed in Fig. [T] 
that the DOS spectrum reproduces energy levels with any symmetry at the quantum level. 
Actually, i^Dos and H are both symmetrical for the exchange of particles, but depends on 
less variables. Thus, the DOS method does not have enough degrees of freedom to study 
quantum states with a particular symmetries and the mass degeneracies due to quantum 
symmetries are not lifted. 

The DOS method presented here shares some similarities with the envelope theory (or 



equivalently the auxiliary field method) [16|, l2ll \24l |25| . For vanishing radial excitations, 



both methods reduce to a similar set of equations: one to fix the energy as a function of a 
mean radius and the other to determine this radius from a transcendental equation. At the 
origin, the envelope theory was developed for D = 3. But, as it relies on the exact solutions 
for harmonic oscillator Hamiltonians, it is not difficult to extend it to any D. Within this 
ramework, the prediction of the envelope theory for the eigenvalues of Hamiltonian ( 129|) is 



2l| 



E 2 = 12 c(2N + L + D) = 12 c(Q + D) . (32) 



This formula is similar to fl30|) . It is an upper bound, but it suffers from a strong degeneracy 
in Q = 2 N + L. Formula ( 130]) does not yield bounds on eigenvalues but finally produces 
better results, since it has the nice feature of lifting the harmonic oscillator degeneracy. 

V. CONCLUDING REMARKS 

The dominantly orbital state (DOS) method is a semiclassical technique to compute 



approximate solutions for quantum eigenvalue problems jl|. It is expected to be valid only 
for high values of the orbital angular momentum and small radial excitations, but it can 
give good results for the whole spectra in some particular cases. This method is very simple 
to implement and needs only the solution of a transcendental equation. In favorable cases, 
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FIG. 1. Eigenmasses E in GeV of Hamiltonian ([21^ for D = 3, and for a = 0.2 GeV 2 and 6 = 0, 
as a function of the number of quanta Q, obtained by an accurate computation in a harmonic 
oscillator basis 10]: All values corresponding to completely symmetrical (circle) and completely 
antisymmetrical (square) states; Only the value corresponding to the mixed symmetry (diamond) 
state for Q = 1. The corresponding approximate values (cross) given by formula (|30p are indicated 
with the quantum numbers (L,N). 



an analytical formula can be obtained. 

Developed at the origin for two-body relativistic systems, it is extended here to treat 
two- and three-body Hamiltonians in D > 2 dimensions with arbitrary kinetic energy and 
potential. A lot of accurate techniques exist to solve numerically two-body problems. So, 
the DOS method is specially interesting for cases in which it can produce an analytical 
solution. But, since quite good results can be obtained by the method, they can also be 
used to check a numerical procedure. 



A three-bod y p rob 



em is always difficult to solve, and techniques to implement it can 



be very heavy lfj, [llj. In this case, the purpose of the DOS method is not to compete 



with these very accurate techniques, but instead to provide rapidly approximate, and even 

11 



analytical, results. Despite the strong constraint imposed on the motion of the particles, 
the results obtained seem reliable: We reproduce the harmonic spectrum exactly and the 
Regge trajectories of relativistic systems with a linear confinement. An interesting feature 
of the DOS method is that it lifts the harmonic oscillator degeneracy which is present in 



other methods 



16 



21 



24 



25l | . This property can be particularly interesting when applied 



to hadronic physics for example 26) . Moreover the DOS method can also be used to check 
a - often complicated - numerical procedure. 
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